Sensitivity to a-variation in ultracold atomic-scattering experiments 
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We present numerical calculations for cesium and mercury to estimate the sensitivity of the 
scattering length to the variation of the fine structure constant a. The method used follows ideas 
Chin and Flambaum [Phys. Rev. Lett. 96, 230801 (2006)], where the sensitivity to the variation 
of the electron to proton mass ratio, /3, was considered. We demonstrate that for heavy systems, 
the sensitivity to variation of a is of the same order of magnitude as to variation of p. Near narrow 
Feshbach resonances the enhancement of the sensitivity may exceed nine orders of magnitude. 
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I. INTRODUCTION 



approximation, the phase shift reads [23 
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Theories unifying gravity with other interactions sug- 
gest the possibility of spatial and temporal variation of 
fundamental physical constants, such as the fine struc- 
ture constant, a = /he, and the electron to proton 
mass ratio, /? = rrie/mp [l|. The search for such varia- 
tions has received considerable interest in recent years, 
and is being conducted using a wide variety of meth- 
ods (see, e.g., reviews @, Q). Two major directions of 
search for variation in a and in /3 are observational stud- 
ies, such as analysis of high resolution spectroscopy from 
dark clouds in interstellar space (see, e.g., Ref. QtMI and 
references therein), and laboratory research, including 
frequency comparison of optical and microwave atomic 
clocks over extended periods of time (see, e.g., [gI. ItI. Hol- 

Great advantage may be gained in these searches by 
focusing on systems that have enhanced sensitivity to 
variation of a or /?. For spectroscopy experiments, such 
an enhancement may occur due to a quasi-degeneracy of 
two close l ying levels of different nature, found in some 
atoms 11, 0, [3111, molecules [13, El, and nuclei 
Chin and Flambaum pointed out that a large enhance- 
ment (on the order of 10^ — 10^^) could also be achieved 
in ultracold atomic collision experiments near Feshbach 
resonances [l^]. In their work, the sensitivity of the scat- 
tering length to variations in the electron to proton mass 
ratio was investigated. Here, we extend this work by fur- 
ther considering sensitivity to the fine structure constant. 

The s-wave scattering length associated with the col- 
lision of two neutral atoms may be written as [ll 



</)= / n-^^-2^i\V{r)-V{^)]dr, (2) 
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where a is the mean scattering length and (f) is the scat- 
tering phase shift. The mean scattering length is fur- 
ther given by a = c(2^C'e/h^)^^'^, with c being a con- 
stant (w 0.47799), n the reduced mass, and Cg the van 
der Waals coefficient characterizing the long-range inter- 
atomic interaction. Employing the semi-classical WKB 



where V{r) is the molecular potential at internuclear sep- 
aration r and the limit corresponds to the inner turning 
point for zero kinetic energy of the colliding atoms, i.e. 
V{Ti) = V'(oo) at the onset of the repulsive wall. 

From Eq. ([T|), we see that the scattering length di- 
verges at values of = {n + |)7r, with n being an inte- 
ger. Such divergences may lead to enhanced sensitivity 
of the scattering length a to variation of the fundamental 
constants: a small change in the phase shift could pro- 
duce a significant change in the scattering length. How- 
ever, the main enhancement is due to the large ratio 
(Ai/me)^/2 - 10^ - 10^ where ^ is the reduced mass 
and ?Tie is the electron mass. This large parameter jus- 
tifies the applicability of WKB approximation, Eq. ([2|) , 
to the motion of a heavy atom in the inter-atomic poten- 
tial, which has an electron-volt scale, and also makes the 
phase very large. As a result, a tiny fractional change in 
the phase cf> produces a dramatic change in tan (0 — f ) 
and the scattering length in Eq. Therefore, the ef- 
fects of the variation of the fundamental constants man- 
ifest themselves mainly via a change in the phase 0. 

From the explicit dependence of 0, Eq. ([5]), on the 
reduced mass /i, it follows that a variation in /3 produces 
an associated change in ip according to the simple relation 
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Here we consider the change of the scattering phase 
shift (f) with respect to variation of the fine structure con- 
stant a. The a-dependence arises from within the molec- 
ular potential V{r) due to rclativistic effects, and the ex- 
traction of this dependence therefore requires ab-initio 
relativistic computations of the molecular potential. We 
compute the potential energy curve for the ground 
state of the Cs2 dimer, and use the previously calculated 
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relativistic and non-relativistic ^T,'^ ground state poten- 
tial energy curves for Hg2 published by Pahl et al. f23 . [25j . 
obtained at the coupled cluster level in the complete ba- 
sis set limit. With these potential energy curves we are 
able to obtain reasonable estimates for the sensitivity of 
(f) to a- variation for Cs2 and Hg2. We present our results 
in terms of the factor Ken where 

5a 



0.000 



54> 



/3- 



(4) 



and Kp = —1/2. From this one can derive the change in 
the scattering length (Eq. ([T])) with respect to variation 
of a and /3 j23 |. 
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The above relation relies on two assumptions. First the 
potential should be deep enough or the reduced mass 
large enough (as for Cs2 or Hg2) to support a large num- 
ber of bound vibrational states, N. Then the phase can 
be estimated as </) w Ntt ^ tt [23|- The second assump- 
tion is in neglecting the variation of the mean scattering 
length, a, which is justified by the fact that it does not 
benefit from the enhanced sensitivity of the tangent func- 
tion appearing in Eq. ([1]). We then go on to extend our 
findings to the case of Feshbach resonances, where even 
greater enhancement of sensitivity to variation of fun- 
damental constants can be expected [2^ . The dynamic 
formation of quasi-bound Cs2 molecules near narrow Fes- 
hbach resonances has been observed by Chin et al. [2^ . 

The scattering length has dimension of length; hence 
we have to specify the units we use. It is convenient 
to use atomic units (of energy, length, etc.). In these 
units the dependence on a appears due to the relativis- 
tic corrections (i.e., terms containing 1/c). Note that 
use of different units would not change our conclusions 
since the variation of the ratio of commonly used units 
(e.g., the Bohr radius and meter) is not enhanced, while 
the variation of the ratio of the scattering length to the 
Bohr radius is enhanced many orders of magnitude; for 
a detailed explanation see Ref. [22i] . 



II. METHOD AND RESULTS 

The first system of interest here is the cesium dimcr. A 
fully analytical potential energy curve for the Cs2 ground 
state has recently been given by Coxon and Hajigeor- 



giou |27|. This curve employs several fitting parame- 



ters to accurately reproduce experimental spectroscopic 
data over 99.24% of the well depth. Additionally, this 
curve has the appropriate asymptotic behavior describ- 
ing the neutral atoms interacting at long-range, namely 
V{r) — >• V{oo) — Ce/r^ as r — > oo. This analytical po- 
tential energy curve, denoted as MLR3 by the authors, is 
displayed in Fig.[TJ Performing the integration in Eq. ([2|), 
we find the scattering phase shift 
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FIG. 1. (color online) Potential energy curve for the 
ground state of Cs2. The circles correspond to values calcu- 
lated as described in the text and the solid line corresponds 
to the VMLRsir) potential of Ref. [13]. The dashed line repre- 
sents the asymptotic part of the true potential energy curve, 
— Ca/r^. All units are atomic units. 



We have computed the molecular potential curve at in- 
cremental values of r using the relativistic molecular four- 
component program package DIRAC [2^ . Details of this 
calculation are saved for the Appendix. Our resulting val- 
ues are overlaid on the MLR3 potential curve in Fig. [1] 
From inspection of this figure, we note a good agreement 
between our computed values and Vmlr3 (f) in the vicin- 
ity of the equilibrium distance re as well as at shorter 
distances. From the equilibrium distance we continue to 
have reasonable agreement up to r w 1.5re « 14 a.u., 
beyond which we compute values quite different to the 
asymptotic behavior of Vmlr3('')- Moreover, the well 
depth is 10% too large. This is due to limitations in the 
basis set applied and the neglect of higher order contri- 
butions in our coupled-cluster expansion. Further, the 
Breit interaction may also change the potential energy 
curve in this region. 

In principle, we may calculate </> using the potential 
energy curve determined from a fit of our computed data 
points, VccsD(f)- We find, however, that this gives an 
entirely unsatisfactory result due to the inaccurate long- 
range behavior of this curve. Firstly, the excessive well- 
depth implies that integration over r in even the short- 
range is erroneous (the integrand of Eq. ^ contains the 
difference V{r) — V{oo)). Secondly, the wrong asymp- 
totic behavior of VccsD('r) implies that a considerable 
amount of phase is inappropriately accumulated in the 
long-range beyond r ?s 14 a.u. (the slow convergence is 
effectively magnified by the square root in the integrand 
of Eq. Clearly, to obtain a reasonable value of (f) it 

is necessary to satisfactorily treat the long-range part of 
the potential. We choose to do so by replacing the er- 
roneous long-range part of our computed curve with the 
proper asymptotic behavior, — Cg/r^, at a suitably cho- 
sen switching point r^. There is a certain degree of arbi- 
trariness in this choice; we set = 14, as at this point 
the deviation of the MLR3 curve from the asymptotic 
curve, —Cq/t^, is roughly equal to its deviation from the 
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computed energy curve, but of an opposite sign, leading 
to error cancellation. Taking the experimental potential 
depth and Cg parameters [27| , we compute the scattering 
phase shift to be 



477, 



(7) 



this result being 2% lower than the value obtained from 
the accurate MLR3 curve, Eq. ([6|). 

We recompute the potential energy curve for differ- 
ent values of a in the neighborhood of ao = 1/137.036, 
and obtain the scattering phase using the method out- 
lined above. It should be noted that the tail part of 
the potential has a considerable contribution to the cal- 
culated phase (for the MRL3 curve, as much as 28% 
of the phase accumulates after r = 14 a.u.). Thus it 
is important to also take into account the variation of 
the Cg parameter with respect to variations in a. We 
obtain its dependence on a in the Appendix, and ad- 
just the integration parameters accordingly. The poten- 
tial well depth, Dg, is rescaled using the following re- 
lation: £ie(a) = [d"^'?'^- lD^^^^{aQ)\D^^^^{a), where 
D^^^°(a) is the calculated depth of the potential en- 
ergy well for a given value of a. By taking the numerical 
derivative of with respect to a we are then able to 
obtain our estimated factor Ka in Eq. (|4]). 

A different strategy was adopted for the mercury 
dimer. High accuracy calculations of both relativistic 
and non-rclativistic potential energy curves were recently 
carried out with the correct long-range behavior from 
symmetry-adapted perturbation theory [H, [2^ . We in- 
tegrate over these potentials to obtain 4ia=ao 4>a=o 
for relativistic and non-relativistic energy curves, respec- 
tively. Assuming linear dependence of molecular proper- 
ties on a^, we use the two points to derive for Hg2. 
Note that we verified that the assumption of the linear 
dependence on works very well for Cs2. In mercury 
the non-linear corrections will be larger; however, they 
should not significantly exceed the natural relativistic pa- 
rameter Z'^o? ^ Z being the atomic number. 

The resulting values for the two systems are 



A'c.(Cs) ^ 0.04, 



^«(Hg) 



-0.3. 



A typical value of relativistic corrections is Z'^o?^ 
which is equal to 0.34 for Hg. This naturally explains 
the magnitude of for this system. For Cs the result 
is smaller, due to an accidental cancellation between the 
contributions of the main and the asymptotic parts of the 
potential energy curve. It should be stressed here that 
the number of approximations employed in this work do 
not allow us to set reliable error limits on the obtained 
values of Ka- However, performing test calculations with 
reasonable variation of the free parameters (e.g., Va and 
Cg) does not change the final result by more than 30%, 
which gives us a rough error estimate. We can never- 
theless conclude that for heavy relativistic systems the 



sensitivity of the scattering phase to a is of the same or- 
der of magnitude as its sensitivity to the electron to pro- 
ton mass ratio [K^ = —0.5). This means that we have 
strongly enhanced sensitivity of the scattering length to 
the variation of a and may use the equations presented 
in the Introduction to estimate this sensitivity. 

A general estimate of the scattering length sensitivity 
for arbitrary atoms is 
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The enhancement is due to the large number N of the 
vibrational levels in the interatomic potential. For Cs, 
where N = 150, a = 280, and a = 95 [13, the enhance- 
ment coefficient is given by 

^ ^ 800(0.5^ - 0.04^). 

a pa 

More accurate result requires knowledge of the experi- 
mental value of the scattering length. Unfortunately, the 
scattering length for Hg2 is not available from experi- 
ment. 



III. EXTENSION TO FESHBACH 
RESONANCES 

In real scattering systems, the molecular potential has 
multiple channels. When a bound state in a closed chan- 
nel near the dissociation limit is tuned close to the scat- 
tering state in an open channel, Fcshbach coupling be- 
tween the two channels can lead to a resonant enhance- 
ment of the scattering length (following Ref. (22j): 



A = a[l + 
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En — E„ 



(8) 



where A is the scattering length near the resonance, a is 
the off-resonant scattering length, IS.E characterizes the 
Feshbach coupling strength, Eo is the energy of the sys- 
tem (in the case of cold scattering we can put £^o = 0), 
and Em is the energy of the bound state. The main effect 
comes from the variation of the small energy denomina- 
tor. Variation of Eq. ^ is then ^ 
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where the small effect of the variation of AE is neglected. 

To explore the dependence of ^ on a and /3, we express 
SA/A in a manner similar to Eq. (|4]), 
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Chin and Flambaum j22| have derived an expression 
of the sensitivity of A with respect to the electron to 
proton mass ratio by exploring the dependence of the 
bound state energy Em on /? and by using the following 
relation between the phase, (pM, and Em, derived from 
the WKB approximation similar to Eq. ([2]), 
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h-^^2ti[Em-Vc.{r)]dr = Mn. 



(10) 



Here, M is the vibrational quantum number of the bound 
state, and and are the classical inner and outer 
turning points, respectively, for a total energy Em- They 
further obtained 12211 . 
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where p(Em) is the density of states at the energy Em 
and can be estimated from the energy splitting, D, be- 
tween adjacent vibrational levels, p{Em) ~ l/D. From 
similarity of expressions (jlOp and ^ , we may provide a 



rough estimate, Ra/Rp ^ Ka/Kp ^ Z a . This means 
that for the heavy atoms the enhancement coefficients for 
the effects of /i and a variations are comparable. 

According to Ref. [22|, this enhancement near nar- 
row Feshbach resonances (e.g., near ij-wave resonance 
in CS2) may reach 10^ — 10^^. Recently, a method for 
measurement of the scattering length to the accuracy 
10~^ was proposed [2^. A combination of these two 
ideas makes the repeating measurements of the scattering 
length a promising method to study the variation of the 
fundamental constants in the laboratory. Here we men- 
tion that the creation of ultracold ^^'^Cs2 from Feshbach 
resonances from Bose-Einstein condensates was demon- 
strated already by Herbig et al. [sO] • In principle, cosmic 
data about scattering lengths and chemical reaction rates 
in the early Universe may be used too. 



IV. APPENDIX: COMPUTATIONAL DETAILS 

The potential energy curves of Cs2 were calculated us- 
ing the 2010 version of the DIRAC program [2^. In order 
to reduce the computational effort, the four-component 
Dirac-Coulomb Hamiltonian was replaced by the infinite 
order two-component relativistic Hamiltonian obtained 
after the Barysz-Sadlej-Snijders (BSS) transformation of 
the Dirac Hamiltonian in a finite basis set [3l|. This 
approximation includes both scalar and spin-orbit rela- 
tivistic effects to infinite order, and is one of the most 
computer time efficient and accurate approximations to 
the four-component Dirac-Coulomb Hamiltonian. 

For Cs2 electron correlation was taken into account us- 
ing closed-shell single-reference coupled-cluster theory in- 
cluding single and double excitations (CCSD). The Faegri 
dual basis set [s^l was used, consisting of 23sl9pl3(i4/2(; 
Gaussian orbitals. Virtual orbitals with energies above 
45 a.u. were omitted, and the 54 outer-most electrons in 



the valence space were correlated. The molecular energy 
curves were recomputed for different values of a, and, 
together with the computed sensitivity of and Cg to 
a, are used to obtain the dependence of the scattering 
phase on the fine structure constant. 

In order to obtain the dependence of the Cq coefficients 
on a we use the following scheme. The van der Waals Cg 
coefficient may be written in atomic units as [33l - [35j 



duj[ap{iijj)Y 
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where ap(w) is the dynamic scalar polarizability. For 
atomic state n, the polarizability is given by [36| 
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with D the electric dipole operator and A„/„ = En' — En 
being the energy difference between states. If a single 
intermediate state n' dominates the summation, we may 
then estimate Cg by 
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where we use /q°°(1 + x'^)~'^dx = 7r/4. 

For the Cs ground state, n = 6si/2, the majority 
(96%) of the of the total polarizability is due to the 
n' = 6pi/2,3/2 intermediate states [36j . The fractional 
change in the Cg coefficient can then be estimated using 
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where the first term is the fractional variation in the 
static polarizability, and the second term represents the 
fractional variation in the first allowed transition energy. 
It should be noted that for Cs, we take into account 
both the 6pi/2,3/2 states, using the weighted average 
A = (Ai/2 + 2A3/2)/3 of the excitation energies. 

The static polarizability and the low excited energy 
states were calculated using the Fock space coupled clus- 
ter approach (FSCC) . Due to the high sensitivity of these 
properties to diffuse basis functions, a slightly larger basis 
was used, consisting of 2hs21plAdQfig orbitals. The po- 
larizability was obtained using the finite field approach, 
similar to Ref. [l^l . The calculations were performed for 
different values of a, in the vicinity of a^, and their de- 
pendence on the fine structure constant was determined 
using numerical differentiation. Our calculated static po- 
larizability for Cs with 396.5 a.u. is in good agreement 
with the most precise values from two other groups (399.9 
a.u. [3^ and 396.0 a.u. [S^) and with experimental mea- 
surements (401.0±0.6 a.u. [H). Using Eq. (Hg, we ob- 
tain 
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